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We discuss the prospects to measure galaxy-cluster properties via weak lensing of 21-cm fluctu- 
ations from the dark ages and the epoch of reionization (EoR). We choose as a figure of merit the 
smallest cluster mass detectable through such measurements. We construct the minimum-variance 
quadratic estimator for the cluster mass based on lensing of 21-cm fluctuations at multiple redshifts. 
We discuss the tradeoff between frequency bandwidth, angular resolution, and number of redshift 
shells available for a fixed noise level for the radio detectors. Observations of lensing of the 21-cm 
background from the dark ages will be capable of detecting M > 10 12 h _1 MQ mass halos, but will 
require futuristic experiments to overcome the contaminating sources. Next-generation radio mea- 
surements of 21-cm fluctuations from the EoR will, however, have the sensitivity to detect galaxy 
clusters with halo masses given enough observation time (for the relevant sky 

patch) and collecting area to maximize their resolution capabilities. 



I. INTRODUCTION 

The hyperfine transition of neutral hydrogen at 21 
cm provides a unique source of cosmological informa- 
tion from the epoch of reionization (EoR) and dark ages 
[HE]. It is the target of several ongoing and near- future 
ground-based experiments [3J, as well as more distant 
prospects such as a Lunar-based observatory [3]. Dur- 
ing the two relevant cosmological epochs, the dark ages 
and the EoR, the hyperfine transition is observed in ab- 
sorption or emission, respectively, against the cosmic mi- 
crowave background (CMB). By measuring this signal at 
different frequencies, the large redshift volume of these 
two epochs can be used to generate independent images 
of the spatial distribution of neutral hydrogen at different 
redshifts [5H7]. 

The image of the 21-cm signal from the dark ages 
and/or EoR should appear on the sky as a random 
field described by some statistically isotropic two-point 
correlation function. If, however, that image is dis- 
torted by weak gravitational lensing from foreground 
matter — either the large-scale inhomogeneous distribu- 
tion of mass in the Universe or by discrete objects, like 
galaxy clusters — then there may be local departures from 
statistical isotropy induced [5HH] . Measurement of these 
local departures from statistical isotropy may thus al- 
low a measurement of the distribution of this intervening 
matter. 

The effects of gravitational lensing of 21-cm 
anisotropies are the same as the analogous effects 
on CMB fluctuations and can therefore be analyzed 
with the tools developed for lensing of the CMB [T5TH7] . 
However, the extension of 21-cm fluctuations to far 
smaller angular scales (limited in principle only by the 
baryonic Jeans mass [5]) than CMB fluctuations (which 
are suppressed on small scales by Silk damping 18 ), 
as well as the possibility to see images of the 21-cm 
background at multiple redshifts, makes 21-cm lensing 



far more promising, ultimately, for weak-lensing studies. 

High-resolution imaging of the mass distribution in 
galaxy clusters using 21-cm lensing was investigated in 
Ref. [10] (and with simulations in Ref. [H]). This work 
concluded that while forthcoming experiments may have 
the potential to provide some initial detections of galaxy 
clusters, the full promise of the technique will likely have 
to await subsequent generations. 

The goal of this paper is to re- visit lensing of the 21- 
cm background by galaxy clusters with an analytic treat- 
ment aimed primarily to help understand the dependence 
of the dctcctability of the signal on experimental parame- 
ters. The aim will be to clarify the experimental require- 
ments for such detections and to assist in the design of ex- 
periments to make such measurements. More specifically, 
we use as a figure of merit the smallest galaxy-cluster 
mass detectable by a given experimental configuration 
and then investigate the dependence of this threshold 
mass on the experimental configuration. 

The plan of the paper is as follows: In Section II, we 
review the 21-cm signal and examine its dependence on 
the observation frequency and on the bandwidth over 
which the signal is observed. In Section III we discuss 
the noise power spectrum of radio interferometers and 
study its dependence on frequency, bandwidth, angular 
resolution, observation time, and collecting area. In Sec- 
tion IV we review how lensing of the 21-cm background 
by galaxy-cluster masses is accomplished. We present a 
quadratic estimator for the weak-lensing convergence and 
derive the noise with which it can be measured. In Sec- 
tion V we construct a minimum-variance estimator for 
the galaxy-cluster mass obtained from the lensing con- 
vergence and derive a formula for the smallest detectable 
galaxy-cluster mass as a function of the various experi- 
mental parameters, assuming an NFW mass profile. In 
Section VI we show the results for galaxy-cluster mass de- 
tectability with next generation interferometers as well as 
futuristic ideal experiments and discuss different trade- 
offs between the experimental parameters. In Section VII 
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we compare the prospects for mass measurements from 
lensing of 21-cm fluctuations with other mass measure- 
ments, mention possible improvements to our estimator 
and discuss some relevant subtleties. We conclude in Sec- 
tion VIII. 



II. 21-CM EMISSION/ ABSORPTION SIGNAL 

We begin by reviewing the physics responsible for pro- 
ducing the 21-cm signal from the dark ages and the EoR, 
whose relative comoving volumes are illustrated in Fig. [T] 
Defining the spin temperature T s as the excitation tem- 
perature of the hyperfine transition (characterizing the 
ratio between the number densities of hydrogen atoms 
in the excited and ground-state levels), the rest-frame 
brightness temperature of a patch of the sky is given by 
Ti, = 7cMBe~ r + Tg(l — e~ r ), where the optical depth 
for the hyperfine transition is |19j 



3c 3 /L4io nm 
I6ku$ T s H{z) ' 



(1) 



A 10 = 2.85 x 10 5 s _1 is the Einstein coefficient for the 
transition, vq — 1420 MHz its rest-frame frequency and 
riui is the local neutral hydrogen density. The brightness 
temperature in this patch, at an observed frequency v 
corresponding to a redshift 1 + z = Vq/v, and the CMB 
is given by the difference 



Ts — Tomb 
1 + z 



(2) 



Hence if the excitation temperature T s in a region differs 
from that of the CMB, the region will appear in emission 
(if T s > Tcmb) or absorption (if T s < T CMB ) against the 
CMB. 

During the dark ages, there is a redshift period 30 < 
z < 200 where neutral hydrogen should be visible in ab- 
sorption against the CMB, as the spin temperature is 
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FIG. 1: An illustration of the comoving volume of the uni- 
verse demonstrating the huge potential for extracting valu- 
able information from 21-cm radiation in the dark ages and 
the reionization epochs. 



coupled to the gas temperature (and cools adiabatically 
as (1 + z) 2 ) while collisions are efficient and drops below 
the CMB temperature (which only cools as (1 + z)), a 
process that peaks at z ~ 70 and lasts until the Hubble 
expansion renders collisions inefficient and T s ~ Tcmb 
again at z ~ 30. The signal from this early epoch is 
not affected by nonlinear density structures nor contam- 
inated by astrophysical sources, which have yet to form. 

A 21-cm signal is also accessible during the EoR, 
at redshifts 7 < z < 13, when newly formed struc- 
ture heated up neutral hydrogen but before the hydro- 
gen became fully ionized. Predicting the brightness- 
temperature signal in this epoch is much harder, as as- 
trophysical noise sources are substantial and our uncer- 
tainties as to the beginning and duration of this period 
are significant. Unlike the dark ages signal, in order to 
plot the signal during reionization we need to trace the 
redshift behavior of the neutral gas fraction which deter- 
mines the optical depth, Eq. |l]), for the hyperfine tran- 
sition. 

In Appendix A, we review the power spectrum of in- 
tensity fluctuations of 21-cm radiation. Eqs. (A3)-(A5) 



describe the angular power spectrum and its approxi- 
mated form in the limit of large and small scales. In 
observations of these fluctuations, an important factor 
is the damping that results from line-of-sight averaging 
in a width Sr (corresponding to an observed bandwidth 
Ais) around the distance r to the desired frequency is, 
whose scale is determined for a given radio interferome- 
ter through the relation [5j 16] 



Sr/r ~ 0.5(Ai//i/)(1 



-1/2 



(3) 



In Fig. [2] we plot the 21-cm dark ages signal power spec- 
trum, Eq. ( |A3[ ), in several redshifts using CAMB [7] in 
the limit of sharp frequency bandwidth and with a band- 
width of Av = 0.1 MHz (resulting in damping above 
I ~ 10 3 ) in the linear regime. We see that for a given 
bandwidth, the signal is roughly within the same order 
of magnitude for most redshifts up to very small angular 
scales, corresponding to arcsecond resolutions. 

In Fig. [3] we plot the 21-cm EoR signal power spectrum, 
assuming reionization is instantaneous and complete at 
z=7, using the approximations in Eqs. (A4) and (A5| 
for a bandwidth of 1 MHz. The damping scale is again 
approximately I ~ 10 3 (as the redshift of reionization is 
an order of magnitude smaller) , and the signal is roughly 
constant for small scales up to > 10 4 , corresponding to 
slightly better than arcminute resolutions. 

In the next section we compare the 21-cm signal with 
the noise power spectrum estimation of radio interferom- 
eters, before heading on to discuss the reconstruction of 
galaxy clusters from the weak lensing signal of the 21-cm 
radiation. 
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FIG. 2: 21-cm power spectrum for different redshifts during 
the dark ages, calculated using CAMB Sources [7 . Solid lines 
are calculated in the narrow-bandwidth limit and the signal 
is damped by the effect of baryon pressure at / > 10 6 . Dashed 
lines are for Av = 0.1 MHz, where the signal is damped due 
to line-of-sight averaging over the bandwidth above Z > 10 3 
as predicted by Eq. Q for these redshift ranges. 
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FIG. 3: 21-cm EoR power spectrum at redshift z — 7, as- 
suming co mple te an d ins tantaneous reionization before 2 = 7, 
using Eqs. ( A4 1 and ( A5 1 (with an amplitude 5Tb{v) ~ 20 mK 
[5]). The solid line is the signal power spectrum for a band- 
width of Av = 1 MHz, which is damped around I < 10 3 
for this redshift according to Eq. j3). Dashed lines are noise 
power spectra for SKA (light gray), SKA with lOxto (medium 
gray) and SKA with 10 x to and 4 x / COVC r (dark gray). We 
see that with the final setup we can reach the maximum 
/ observable with SKA's baseline of 6 km which is roughly 

fcover ^ 2 X 10 . 



III. EXPERIMENTAL PROSPECTS 

Following Ref. [5], we examine the noise power spec- 
trum of a radio interferometer array. We denote by 
Z cover (z/) = 2ttD/X the maximum mode at frequency 
v (corresponding to wavelength A) that can be mea- 
sured with an array of dishes with maximum baseline 
D covering a total area A to tai with a covering fraction 
/cover = NdishAlishMtotal- For uniform Fourier area cov- 
erage, a system temperature T sys , a frequency window 



A;/, and an observing time t , the noise power spectrum 
is 



l 2 C? 



(27r) 3 r s 2 y » 



2 

cover 



(4) 



Thus, for uniform Fourier coverage in an experiment tar- 
geting a frequency v with bandwidth Av, we conclude 
that C? - const = 27r/3(^)(/ c 2 ovcr toAv)' 1 , where fi{v) = 
(2i:) 2 T 2 ys (v)[l covcr (v) 2 ]' 1 . As we can see from Figs. [2] 
and |3j on the relevant scales for cluster reconstruction, 
the 21-cm radiation power spectrum at a given band- 
width behaves roughly as l 2 Ci ~ const = 2na(v, Av). 
To find the value ? max beyond which the noise Cf is no 
longer negligible, we compare C/ max = C" to find 



I 2 = 

max 



(/cover *0Al/). 



(5) 



Observing the 21-cm anisotropies during the dark ages 
will be very challenging. The redshift range 30 < z < 200 
of absorption of the CMB at 21-cm {y = 1420 MHz) in 
the rest frame of the neutral-hydrogen gas corresponds 
to very low frequencies 7 MHz < v < 46 MHz, where 
the sky temperature (dominated by foreground sources 
of synchrotron emission from the galaxy and from ex- 
tragalactic sources) in regions of minimum emission at 
high galactic latitudes, approximately given by T sys ~ 

180 (zVl80MHzr 2 - 6 K : reaches - 10 4 - 10 6 K, many or- 
ders of magnitude above the signal ~ 1 mK (as seen in 
Figs. |] and [3}. 

While both galactic and extragalactic emissions vary 
smoothly and could be possibly subtracted by taking 
observations at two closely spaced frequencies, the ad- 
ditional sources of interference are harder to overcome. 
Terrestrial radio frequency interference is abundant in 
this range and is not spectrally smooth. Another source 
of contamination is the ionosphere which causes phase 
distortions in the cosmic signal and turns opaque at fre- 
quencies v < 20 MHz (corresponding to z > 70). A futur- 
istic experiment, based on placing a dark ages observa- 
tory on the far side of the Moon (which has no permanent 
ionosphere and its far side is shielded from terrestrial ra- 
dio interference) was suggested in Ref. [4] . With a base- 
line on the order of > 10 — 100 km, angular resolutions 
corresponding to Z max ~ 10 4 — 10 5 (for source redshifts 
z = 30 — 300) might be reached in such an experiment, 
compensating for covering fraction and system temper- 
ature with prolonged observation time. With sufficient 
frequency coverage, the huge volume of CMB absorption 
in neutral hydrogen during the dark ages, 30 < z < 200, 
can be used to beat down the noise (by combining red- 
shift slices, as we describe in the next Section) and mea- 
sure the mass profile parameters of galaxy clusters to 
high accuracy. 

To reach a measurement within the decade, the best 
candidate is the currently planned next generation exper- 
iment to measure 21-cm emission from neutral hydrogen 
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during the EoR, the Square Kilometer Array 1 (SKA). 
With a design based on an extended region of Z?~6 km 
where Z covor (z/) ~ 10 4 for the relevant frequencies in the 
epoch of reionization and / CO vcr ~ 0.02, SKA can reach 
an observing time of order to ~ 1000 hrs in one season 
covering up to an area of ~ 2?r sr in the sky. 

In Fig. [3]we plot several noise power spectra matching 
the current plans for a three-season run of SKA measur- 
ing at a frequency corresponding to redshift z = 7 with a 
bandwidth of 1 MHz, as well as scenarios with an order of 
magnitude larger observation time, and with four times 
the coverage fraction. With the current plan, the maxi- 
mum scale observable with a signal to noise greater than 
one is only 0(l max ~ 10 3 ). To reach l max > 10 4 will re- 
quire a considerable increase in coverage fraction and / or 
observation time. 



spectrum Ci defined by 



(8) 



The Fourier coefficients for the observed intensity are re- 
lated to those of the unlensed intensity and deflection 
field by 



V 



(9) 



from which it follows that for a fixed deflection field, the 
observed intensity satisfies 



h k 



I?Ir) = SlCiSr 



h-k 



> E L ■ (hC h + l 2 C h 



(10) 



IV. LENSING OF 21-CM FLUCTUATIONS BY 
A GALAXY CLUSTER 



where L = l\ + 1% 
Since 



-L 2 4> z /2, each l x -l 2 pair of measured in- 



tensity Fourier modes with wavevectors l\ + I2 = L pro- 
vides an estimator, 



We now review the distortion to the 21-cm fluctuations 
induced by weak gravitational lensing. We assume that 
the lensing distortion takes place over a relatively small 
region of the sky, as should occur for lensing by a galaxy 
cluster, so that we can work in the flat-sky limit, where 
the analytic expressions are simpler. 

Let Iq(9) be the 21-cm intensity at position 9 on the 
sky. Lensing will deflect photons from by an amount 
89(9) = V 4>{6). Here, (f)(9) is the projected potential, re- 
lated to the convergence k(9) = 2£(0) /S cr , by V 2 ^ = 2k, 
where is the surface mass density of the interven- 
ing cluster, and S cr = AirGDdDd s /(c 2 D s ) is the critical 
surface mass density in terms of the observer-lens, lens- 
source, and observer-source angular-diameter distances 
Dd, Dds, and D s , respectively. The observed intensity is 
thus 1(9) =I (9 + 69) ~ I Q (9) + 89- VI (9). 

We suppose that intensity is measured over some 
square patch of sky of solid angle fi, decomposed into 
jV pix pixels at positions 9 i7 for i — 1, 2, . . . , JV p jx, sur- 
rounding the cluster. The intensity can then be written 
in terms of Fourier coefficients, 



Ir = 



n 



(6) 



L 



v- 



2 L ■ (hC h + l 2 C h ) 
and the variance of this estimator is 

l 4 c™ ap c™ ap n 2 




L-(hC h +l 2 C h ) 



(11) 



(12) 



where C'[ nap = C; + Cf is the power spectrum of the 
intensity map, including the noise Cf. 

We can then sum the estimators n l i' 12 over all h + h = 
L pairs (correcting for double counting of triangles with 
I L — I) with inverse-variance weighting to obtain the 
minimum- variance estimator, 



E 



L 2 



? L-{lC l + (L-l)C^_ A ) 



[L ■ (Id + (L- l)C l£ _ n 



-M2 



(13) 



and 



1(9,) = 



r 



where the sum is over the N p { x Fourier modes I. 



(7) 



The 



two-point intensity correlations are described by a power 



1 http://www.skatelescope.org 



NZ- 



2n s -^[L-(ic l + (L-i)c ]l _ A )] 2 



T 



(14) 



is the noise power spectrum for Kg. Equivalcntly, 
(|«-J| 2 ) = (2n) 2 S(0)N £ = nJVTg is the variance with which 
can be measured. 
We now use the continuum limit il -1 Y]f — > 
J d 2 l 7(27r) 2 and employ the simplifying assumption, 
Eq. ([5]), that the noise is small, Cf <C Cj, so that C ; map ~ 
C/, up to some scale i max - Then, in the limit L <C I we 
approximate \L — l\ ~ I — Lcoscj), where cos0 = L ■ I, 
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and then to first order, C|£_n ~ Ci — L(cos(/))(dCi/dl), 
which yields 



decreased, which limits the actual number of slices that 



can be combined in Eq. (17). 



NZ 1 



L 4 



I 2 

max 
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l 



L • ICi + L • (L - l)C, 



\l-i\ 



2Ci C|f_f] 
91nC, 3 /<91nC, 



<91nZ 8 V dhil 



(15) 



Using the approximation, Eq. the noise for a single 
slice in an experiment at frequency v with bandwidth Av 1 
a maximum baseline corresponding to ^ C over, a coverage 
fraction f covei , and observation time to, is approximately, 



47T 

P 
max 



4tt/3(i/) 
a(^,A^)(/ c 2 ovcr t Ai.)- 



(16) 



We can increase our signal considerably by changing 
the frequency at which the 21-cm map is made and 
thereby focus on spherical shells of neutral hydrogen at 
different redshifts. The above results for a single redshift 
slice can be extended to make use of the full redshift vol- 
ume [5] by discretizing the z direction into radial compo- 
nents so that C^k is the power in a mode with angular 
component I and radial components k = 27rj/i?, where 
R is the total radial length of the volume. Under the 
assumption that modes with different k are independent, 
a total-volume estimator is built by summing the indi- 
vidual estimators, and the corresponding noise variance 
is 



L 4 /4 



E 



21 



L- (L- 0C|L-i],/J 2 



n^map^map 
ZU l,k °|f-il,fc 



JL • lG(u) + h .(L -l)C [£ _ n {v)f 



(17) 



where we switched to a sum over the corresponding fre- 
quencies in the second line. In general, to estimate the 
number of slices available [Bj, we can think of 21-cm maps 
at two different frequencies that correspond to spherical 
shells separated along the line of sight by a comoving 
distance Ar. These will be statistically independent at 
the highest I provided that (Ar/r„) > l^ 1 . An experi- 
ment that covers a spatial range 5r or a frequency range 
v\ — vi around a frequency v will yield a total number of 
N z ~ (Sr/Ar) ~ l(Sr/r u ) ~ 0.5i((i/i - ^)/f)(l + z)- 1 ' 2 
statistically independent maps, which is roughly N z ~ 
1500 for the frequencies corresponding to the EoR at the 
maximum resolution of SKA and up to N z ~ 5000 for 
a dark ages observatory with a baseline of 100km (here 
we neglect the fact that different frequency bins are also 
somewhat correlated by contaminations). In practice, 
however, the tradeoff is a complicated one. While the 
signal is larger for narrow bandwidths, we saw in the 
last Section that the noise power spectrum around a cer- 
tain frequency increases as the frequency bandwidth is 



ESTIMATING CLUSTER PROPERTIES 



Applying the estimator in Eq. ( 14 1 for the convergence 



/?£ to a patch of sky around a galaxy cluster (assuming 
we know the location of its center to enough accuracy, say 
from other sources such as SZ surveys), we can retrieve a 
2D image of the weak-lcnsing convergence of the cluster 
which we can study versus theory or simulation. The to- 
tal number of modes is set by the total number of pixels, 
and we denote each mode by Li. If Tw~. is the measured 
value for the pixel i with variance (|kJ~| 2 ) = and 
the corresponding theoretical value is k^ 1 (Msd) (calcu- 
lated for some fiducial value A/fid), then an estimator for 
the halo mass per pixel is given by 



Mi 



KfiM Rd y 



M, 



fid ■ 



with variance 



(18) 



(19) 



From these we can build a minimum-variance estimator, 
EMJdMl 2 ) 



M 



EV<IW> 



(20) 



for the mass as a weighted sum over the pixels. The 
variance of this estimator is, under the null hypothesis, 



Mi d QN £ 



N I 2 

Mf d fi47T 



N I 2 

1 v 2 'max 



5>| h (M fid )| 2 |W £ | 



'ficT 
N I 2 
M 2 d 4n 



d 2 L 



Kf(M 6i )W £ \ 



M 2 d 4TT J (2tt) 2 ' l 



d 2 



d 2 (pW{<p)K{\e - <p\) 



(21) 



where in the second line we substituted 



L N I 2 



3 L 2 /2Lf llax _ 



47T 



NP (22) 



with the exponential included to describe roughly the 
transition between those L modes that can be measured 



and those that cannot. In the third line in Eq. (21) we 
took the continuum limit J^l ^ ^ / 'PL/Qii) 2 , an d in 



G 



the fourth we used Parseval's theorem to switch to real 
space, resulting in a convolution of the convergence with 
a two-dimensional Gaussian filter, 



W{<p) = 



1 



27T<9? 



(23) 



with smoothing scale 9 S given by 6 S = 7r/Iy max (we as- 
sume here that we can push L max close to Z max which is 
verified numerically). Assuming a spherically symmetric 
profile and cutting off the integral at A when the signal 
becomes negligible, we get 



M 



N P 
2M| d 



ede 



-, 2 



(fdifW(ip) / d(f)K(^0 2 + ip 2 + 26ip cos < 



(24) 

This equation, together with Eq. (5j), allows a straight- 
forward examination of the capabilities of different tele- 
scopes to reconstruct a given lensing source. The focus 
of this work is detection of a galaxy cluster, but a sim- 
ilar formula can be used to estimate the signal-to-noise 
for reconstruction of model parameters for other lensing 
sources, as we discuss in the Conclusion. 

The remaining task is to calculate the convergence pro- 
file of a given cluster. We model the mass profile of the 
galaxy cluster by an NFW profile [2"3"] . 



p(r) 



'r/r s (l + r/r s ) 2 ' 



(25) 



where the scale radius r s and normalization p s are often 
described by the concentration parameter c = r v ; r /r s and 
the cluster mass M = 47rrf p s [ha(l + c) — c/ (1 + c)] . Here , 
r v ; r is the radius within which the enclosed mass M is 200 
times the average mass of the same volume in a critical 
density universe. 

The convergence k(9) of the NFW profile is given by 

K (6 = r/D L )= 2 ^f(0/(r s /D L )), (26) 

2jci- 

where the functional form of the projected mass density 



fix) 




vT 5 



arctan 



arctanh 




In Fig. [4] we plot the convergence of a cluster at red- 
shift z = 1 with mass of Mi = 5 x 10 14 /T^M© and 
concentration parameter c = 3 for a source at redshift 
z = 7. Plugging this in the integral of Eq. ( 24 ) , we find 



that the signal-to-noisc is saturated at the virial radius, 
so that effectively A = r v „ , an order of magnitude above 
the corresponding Einstein radius. We also find that the 
dependence on the concentration parameter is small. 
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FIG. 4: The radial convergence of NFW clusters with mass 
of M — 5 x 10 14 /i _1 Mq and concentration c = 3 located at 
redshift z = 1, where the 21-cm source is at redshift z = 7. 
The vertical line marks the virial radii of the clusters where 
the signal becomes negligible. The two shaded regions show 
the amount of signal integrated with the ideal resolutions of 
SKA at 1 arcmin and a dark ages observatory with a 100 km 
baseline at ~ 6 arcsec. 



VI. RESULTS 

To estimate the signal-to-noise with which this cluster 
mass can be reconstructed using the weak-lensing sig- 
nal of 21-cm radiation measured by a radio interferome- 
ter, we consider the expected signal from the EoR from 
SKA (as described in the previous Section) , an upgraded 
SKA with four times the collecting area, both with a 
bandwidth of 1 MHz, and an ideal experiment (capable 
of sustaining Cf > Cf up to ? max = Z C ovor) with the same 
resolution as SKA (Z COV or ^ 10 4 ) with a bandwidth small 
enough to reach the maximum number of independent 
redshift slices. We also include extremely optimistic lim- 
its for an ideal dark ages observatory with a baseline of 
100 km (for which l covcr ~ 10 5 ). 

In Fig. [5] we plot the smallest detectable mass as a 
function of redshift, for different experimental capabili- 
ties. We see that with the current plans for SKA, the 
weak-lensing reconstruction of clusters as considered in 
Fig. [4] is completely beyond reach, unless about an order 
of magnitude more observation time is dedicated to the 
patch containing the target galaxy-cluster. Future ex- 
periments will narrow this gap and enable the mass mea- 
surement of significantly smaller mass halos using 21-cm 
weak lensing. 

These detection prospects involve intricate tradeoffs 
between the different experimental parameters. For ex- 
ample, if we increase observation time by 2, we can reach 
the same / max with half the frequency bandwidth, which 
will also allow us to use twice the number of z-slices, in- 
creasing the S /N by y/2. Alternatively, we could use this 
to reach a larger Z max (paying a small price for the damp- 
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FIG. 5: The la limit for weak lensing reconstruction for clus- 
ters in different redshifts up to the end of reionization (z ~ 7) 
using SKA (thick solid line), SKA with four times the coverage 
fraction or sixteen times the observing time (thin solid line), 
both with a bandwidth of 1 MHz, and for an ideal experiment 
with SKA resolution (dashed line). We also plot the expected 
limits for an ideal dark ages observatory with a 100 km base- 
line (dotted line). 



ing in the signal power spectrum), which also increases 
the S/N, but this will still be limited by the maximum I 
set by the baseline of the interferometer. As mentioned 
in Section IV, one does not always do better by splitting 
up into smaller bandwidth bins because this increases the 
noise, yielding a smaller i max , which reduces the S/N and 
might also leave galaxy-cluster scales beyond reach. In 
addition, as we also discussed in Section IV, there is a 
lower limit to the frequency bandwidth below which the 
slices become correlated. 

Another experimental issue is that as the observed fre- 
quency is increased in order to use multiple redshift slices, 
a larger baseline is needed to cover the same scales. More- 
over, as can be seen for the dark ages in Fig. [2] the de- 
pendence of the 21-cm signal strength on the redshift is 
not trivial as well. 



VII. DISCUSSION 

We have constructed a minimum-variance estimator 
for the mass of a galaxy cluster using the weak-lensing 
convergence reconstruction of 21-cm fluctuations from 
the dark ages and the EoR. As has been suggested 
for CMB lensing reconstruction, possible improvements 
to the results shown here can stem from using mod- 
ified weak-lensing estimators, enhancing sensitivity to 
small scales as in Ref. 24; or iteratively approaching a 
maximum-likelihood estimator as in Ref. 25 . Another 
option which has been discussed for CMB measurements 
(e.g. Refs. [17l HH [25]) is to stack reconstructed images 
of N different clusters to yield a y/N improvement in the 
signal-to-noise for fitting an overall mass profile (much 
like the use of different 21-cm redshift slices of the same 



cluster in Eq. (17)) 



In comparison to other sources for weak-lensing mea- 
surements, such as galaxy shapes [26l[27], or CMB fluctu- 
ations [T5] (detected recently [25]), the potential in 21-cm 
measurements is far greater. Even deep galaxy surveys 
will be limited to lower redshifts with small sky cover- 
age and will not exceed arcminute resolutions. CMB 
lensing reconstruction at small scales is difficult because 
of Silk damping at these scales and the obstruction by 
the Sunyaev-Zeldovich effect on the CMB photons pass- 
ing through the moving clusters. Even with ideal ex- 
periments at arcminute resolutions, the corresponding 
mini mum galaxy-cluster mass detectable according to 
Eq. ((24) with a CMB experiment is M ~ 10 15 h~ 1 M Q . 
An average mass profile of M ~ 10 14 ft- _1 M Q can be mea- 
sured with reasonable signal to noise by stacking thou- 
sands of galaxies, but the detailed reconstruction of an 
individual cluster is beyond the reach of these alternative 
methods. 

Finally, an issue which we have neglected in this work 
is the influence of non-gaussianities which would have to 
be taken into account to yield accurate predictions. Par- 
ticularly important is the effect of nonlinear structure 
[121 [29] , which is relevant during the epoch of reioniza- 
tion, when certain patches of the IGM become substan- 
tially ionized well before its end. In order to properly 
account for these effects, the details of the reionization 
process will have to be uncovered. However, under the 
assumption that these features appear in higher resolu- 
tions than that of our considered interferometers [5J [TU] , 
there will be no connected four-point function contribu- 
tion to the variance of the quadratic lensing estimator 
used here, and this treatment remains valid. 



VIII. CONCLUSION 

When 21-cm fluctuations become accessible to obser- 
vations at small angular scales, the application of weak- 
lensing reconstruction methods will open the door to un- 
precedented precision measurements of local structure. 
We have found here that galaxy clusters can be de- 
tected by lensing reconstruction with futuristic experi- 
ments measuring CMB absorption by neutral hydrogen 
during the dark ages 30 < z < 200 down to halo masses 
of order M > 10 12 /i _1 Mq. Next-generation interferom- 
eters measuring emission from hot neutral-hydrogen gas 
during the epoch of reionization 7 < z < 13 will be lim- 
ited to M > 10 15 hr 1 Mq with the currently planned 
specifications. If, however, the collecting area or the ob- 
servation time can be increased, they may be able to push 
the limit down to M > 10 13 M Q (in particular, this 
can be reached by increasing the observation time ded- 
icated to the target patch alone by two orders of mag- 
nitude). While these are challenging goals that remain 
unattainable in the near future, the potential achieve- 
ments discussed here provide more motivation to invest 
in alleviating the experimental limitations. 

Meanwhile, we can use the prescription described 
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here to construct estimators for model parameters of 
other structures, from standard isothermal spheres and 
voids to more exotic structures such as cosmic textures 
[501 |3"T] or overdensities created by pre-inflationary par- 
ticles 32, 33J . In future work we shall address the weak- 
lensing detectability, with 21-cm weak lensing, of voids 
or textures that might be responsible for the cold spot in 
WMAP CMB data El [35] . 
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Appendix A: Power spectrum of 21-cm fluctuations 

We review the angular power spectrum of intensity 
fluctuations in the 21-cm signal induced by large-scale 
density inhomogeneities during the dark ages and the 
EoR. Following Refs. [H [5], we start in three dimen- 
sions and first define the dimensionless fractional per- 
turbation 6hi(x) = [STb(x) — ST b]/ 5Tb of the brightness 
temperature (where 5Tb is its mean), whose power spec- 
trum in Fourier space is given by (5 h i O^-i) 5 h i (^-2)) = 
(27r) 3 <5 3 (ki + k2)P(fci, z) with its dimensionless counter- 
part A 2 (k,z) = (k 3 /2n 2 )P(k,z). We move to two di- 
mensions by integrating along the line of sight, 

5 HI (n, v)= Jdr W(r, r„) 5 H i(&, r), (Al) 

with a projection window function W(r, r v ) that peaks at 
the radial distance r v to the desired frequency v within a 



width 5r corresponding to an observational bandwidth 
Aza Expanding in spherical harmonics 5m{n,v) — 
ai m (v)Yi m {h) and using the planar wave expansion 

lm 

elk ' x = Y.i m ^ il 3l{kr)Y? m $)Y lm {h) (where j t (x) is the 
spherical Bessel function of order I), we get 



ai m (v) = 



d 3 fc 
(2tt) 3 



5 H i{k,v)ai(k,v)Y l * m (^) 



ai(k,v) = 5T b (v) I drW(r,r v )ji(kr) 



(A2) 



The brightness temperature power spectrum of 21- 
cm fluctuations, defined as (ojro( I/ )o*' m '( I// )) = 
5w5 mm >Ci{v,v'), is then given by 



Ci(v,v') = 4tt j dk A2 ^ ,z) ai (k,u) ai (k,i/). (A3) 

For a pure power-law power spectrum, A 2 (fc,z) 
A 2 (//rv,z) and for large angular scales I5r/r <C 1 the 
line of sight integration is approximately a delta function 
and we get 



l 2 Ci{v,v) „-2, , a2 



2?r 



<x5Tl(v)A 2 (l/r v ,z). 



(A4) 



For the purposes of reconstructing galaxy-cluster profiles, 
we are mostly interested in the smallest observable scales, 
which satisfy the limit I5r/r 3> 1. Applying the Limber 
approximation [5T] in Fourier space [55] yields 



2^ °c 5T b (v)A (l/r u ,z) — , (A5) 



where we see the suppression of small-angle fluctuations 
resulting from averaging out of modes k > 1/Sr in 
Eq. jm 
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